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EFFECT  OF  DIFFUSION  FACTORS  ON  FLAME  STABILITY 


V.  N.  Iyevlev  and  S.  A.  Gol'denberg 


This  article  considers  the  causes  for  the  devia¬ 
tion  from  the  universal  relationship  of  experimental 
values  of  critical  velocities  of  flame  blow-off  from 
small  stabilizers  [1,  2].  The  experimental  investi¬ 
gation  was  carried  out  with  a  gasoline-air  mixture 
using  conical  and  cylindrical  stabilizers. 


Certain  anomalous  phenomena  were  detected  when  carrying  out  cer¬ 
tain  experimental  investigations  into  flame  stabilization  with  poorly 
streamlined  bodies  and  when  elucidating  the  dependence  of  the  flame 
blow-off  velocity  from  stabilizers  on  the  fuel  mixture  composition. 

For  example,  as  early  as  in  William's  study  [3],  a  disconformity 
of  the  position  of  curves  of  flame  blow-off  velocity  and  the  curve  of 
normal  flame  velocity  was  established  (as  a  function  of  the  fuel  mix¬ 
ture  composition  or  excess  air  ratio),  which  was  especially  evident  in 
the  region  of  maximum  values  of  these  velocities.  The  data  derived  in 
these  studies  were  not  governed  by  the  universal  relationship 

:  Tr*(-f)-  <*> 

This  type  of  unconformity,  which  was  observed  in  the  studies  men¬ 
tioned,  could  not  be  explained  by  the  experimental  errors. 
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Therefore,  the  task  of  this  experimental  investigation  was  to 
elaborate  the  described  phenomenon  on  small  stabilizers,  conical  and 
cylindrical  rods. 

The  experiments  were  carried  out  with  a  gasoline -air  mixture. 

The  stabilizers  were  installed  along  the  flow  axis  of  the  fuel  mixture 
at  a  distance  of  5  nm  from  the  nozzle  of  the  burner,  a  tube  J>2  mm  in 
diameter.  We  used  cylindrical  (rod  d  =  4  mm)  and  conical  (d  *  5  mm) 
stabilizers.  As  a  result  of  the  experiments  we  obtained  values  for 
the  flame  blow-off  velocity  as  a  function  of  the  fuel-mixture  composi¬ 
tion  (excess  air  ratio  a). 

It  wa3  revealed  from  the  experiments  that  the  blow-off  velocity 
w,  as  compared  with  the  curve  of  normal  flame  propagation  velocity  un, 
shifts  toward  the  side  of  rich  mixtures:  the  peak  of  the  blow-off 
curve  lies  at  a  =  0.7,  and  the  peak  of  the  curve  un,  at  a  «  0.9.  When 
using  other  fuels,  the  peak  can  shift  toward  lean  mixtures.  The  lat¬ 
ter  was  observed  In  experiments  with  methane-air  mixtures  [6]:  the 
peak  of  the  blow-off  curve  was  at  a  »  1.1,  and  the  peak  of  the  normal 
flame  velocity  at  a  *  0.95* 

The  shift  of  the  peak  of  the  blow-off  curves  [?,  4]  has  been  ex¬ 
plained  by  the  difference  in  the  numerical  values  of  the  diffusion  co¬ 
efficient  of  fuel  D  and  the  coefficient  of  thermal  diffusivity  of  the 
mixture  a.  The  latter  coefficient  for  air  mixtures  of  hydrocarbons  is 
for  all  practical  purposes  determined  by  the  thermal  diffusivity  of 
air,  where  for  gasoline  D  <  a,  and  for  methane  D  >  a. 

•  t — • 

The  most  direct  method  of  elucidating  the  role  of  molecular  dif¬ 
fusion  In  the  flame  stabilization  process  is,  evidently,  a  study  of  the 
composition  of  combustion  products  and  their  temperature  in  the  circu¬ 
lation  zone  behind  the  stabilizer. 
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Corresponding  experiments  were  carried  out  on  a  conical  stabilizer 
(d  »  7  mm)  and  a  longitudinally  streamlined  rod  (d  •  4  vm). 

The  combustion  products  were  gathered  from  the  core  of  the  circu¬ 
lation  zone  close  to  the  end  of  the  stabilizer  by  means  of  a  thin  med¬ 
ical  needle  with  an  outside  diameter  of  1  mm.  The  needle  was  soldered 
Into  the  air-cooling  tube.  The  temperature  In  the  circulation  zone 
was  determined  by  the  sodium-line  reversal  method. 


TABLE  i 

Excess  Air  Ratio  (Initial  Mixture  and  Beyond 
the  Stabilizer)  and  Composition  of  Combustion 
Products  in  Circulation  Zone 
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We  see  from  Table  1  that  the  excess  air  ratio  calculated  from  the 
analysis  of  the  combustion  products  in  the  circulation  zone  beyond  the 
rod  (d  =  4  mra)  differs  appreciably  from  a  of  the  initial  mixture.  This 
difference  was  not  observed  in  the  case  of  the  conical  stabilizer 
(d  =  5  mm). 

If  the  data  on  the  blow-off  velocity  at  the  rod  (d  =  4  mm)  Is 
plotted  on  a  graph  in  relation  to  the  true  values  of  a  in  the  circula¬ 
tion  zone  (Pig.  1),  the  blow-off  curve  is  shifted  to  the  right  and 
will  occupy  a  position  (dashed  line)  corresponding  to  the  position  of 
the  curve  of  normal  flame  velocity.  In  this  case  the  peaks  of  both 
curves  will  be  in  the  region  a  ■  0. 9  -  0.93.  Similar  data  for  the 
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Fig.  1.  Dependence  of  blow-off 
velocity  w  (m/sec)  on  excess  . 
air  ratio  a  of  the  initial  mix¬ 
ture:  l)  for  rod  d  =  4  mm; 

2)  for  conical  stabilizer  d  » 

=  5  mm;  3)  on  the  excess  air 
ratio  in  the  circulation  zone 
for  rod  d  =  4  ram. 


Fig.  2.  Dependence  of  tempera¬ 
ture  in  circulation  zone  £  (°C) 
on  the  excess  air  ratio  a:  i) 
for  rod  d  ■  4  ram  (the  curve  is 
plotted  in  conformity  with  the 
characteristic  of  the  initial 
mixture);  2)  for  conical  sta¬ 
bilizer  d  ■  5  nim. 
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conical  stabilizer  (d  *  5  mm)  did  not  yield  a  divergence  between  a  of 
the  initial  mixture  and  a  of  the  combustion  products  in  the  circulation 
zone,  and  therefore  shifts  of  the  curves  of  blow-off  velocity  are  not 
revealed. 

In  addition  to  determining  the  composition  of  the  combustion  prod¬ 
ucts  in  the  circulation  zone  beyond  the  rod  (d  =  4  mm)  and  conical 
stabilizer  (d  *  5  Dim),  we  measured  the  temperature. 


TABLE  2 

Results  of  Temperature  Measurements  in 
_ the  Circulation  Zone _ 
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lrfe  see  from  Table  2  that  the  temperature  values  In  the  circula¬ 
tion  zone  calculated  by  analyzing  the  combustion  products  are  close 
to  those  directly  measured. 

Temperature  curves  are  plotted  in  Fig.  2  for  the  rod  d  =  4  mm  and 
conical  stabilizer  d  *  5  ddbu 

From  a  comparison  of  Figs.  1  and  2  we  see  that  the  position  of 
the  peaks  on  the  temperature  curves  corresponds  to  the  curves  of  blow- 
off  velocity.  Evidently,  If  we  assume  for  the  basic  characteristic  the 
true  composition  of  the  combustion  products  a  in  the  circulation  zone 
of  the  rod,  the  temperature  curve  for  it  will  be  shifted  to  the  right 
and  the  peak  of  the  curve  is  at  a  point  corresponding  to  a  **  0.9* 

Thus,  direct  measurements  show  that  in  the  circulation  zone  be¬ 
yond  the  rod  d  *  4  mn,  the  composition  of  the  c cxnbus t ion  products 
changes,  the  nitrogen  content  increases  in  them,  and  owing  to  this  the 
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temperature  drops.  In  this  case  the  composition  of  the  combustion 
preducts  after  total  combustion  and  their  temperature  correspond  to 
the  combustion  products  obtained  on  ignition  of  the  initial  mixture 
a  -  1.23. 

In  Table  2  this  is  illustrated  by  the  magnitudes  of  nx  and  na 
which  are  the  ratio  of  the  number  of  carbon  atoms  to  the  number  of 
nitrogen  atoms  in  the  initial  mixture  (nx)  at  a  «  0.9  and  in  the  cir¬ 
culation  zone  n2.  In  the  case  of  the  conical  stabilizer  d  =  5  mm* 
when  no  shift  of  the  curves  to  the  right  is  observed,  the  ratio 
n2  :  nx  =  1.  In  the  case  of  rod  d  =  4  mm  where  there  is  a  shift  of 
the  blow-off  curves  and  the  composition  of  the  combustion  product 
(danai  =  1* 25)  in  the  circulation  zone  does  not  correspond  to  the  ini¬ 
tial  mixture,  the  ratio  is  n2  :  nx  <  1. 

The  described  phenomena  associated  with  the  effect  of  molecular 
diffusion  under  specific  conditions  of  flame  stabilization  were  noted 
in  works  on  small  stabilizers  [6]. 

A  change  in  the  composition  of  combustion  products  in  the  circu¬ 
lation  zone  beyond  the  small  stabilizer  and  also  a  change  of  their 
temperature  in  comparison  with  the  initial  mixture  leads  to  a  change 
in  the  intensity  of  the  chemical  source  at  the  interface  of  the  fresh 
mixture  flow  and  the  circulating  heated  combustion  products.  Therefore, 
in  these  cases  the  use  of  the  value  of  un  for  calculations  Is  not  sub¬ 
stantiated,  it  should  be  specially  calculated.  However,  the  method 
for  such  a  calculation  has  still  not  been  developed. 

Numerous  calculations  show  that  a  sufficiently  good  generaliza¬ 
tion  of  experimental  data  on  blow-off  from  small  stabilizers  can  be 
obtained  if  we  artificially  combine  with  respect  to  a  the  curves  of 
blow-off  velocity  with  the  curves  of  normal  flame  velocity  or  with 
blow-off  curves  obtained  on  large  stabilizers  (Fig.  1). 
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The  certain  possible  error  In  this  case  Is  associated  with  the 
fact  that  the  blow-off  velocity  curves  should  be  shifted  not  only  with 
respect  to  a  but  also  with  respect  to  temperature.  This  Is  due  to  di¬ 
lution  of  the  mixture  In  the  circulation  zone  as  a  result  of  diffusion 
of  nitrogen  (air)  from  the  Initial  cold  mixture. 

The  authors  express  their  sincere  thanks  to  L.  N.  Khitrin,  cor¬ 
responding  member  of  the  USSR  Academy  of  Sciences  for  his  very  valu¬ 
able  help  in  carrying  out  the  study. 

DESIGNATIONS 

w  Is  the  flow  velocity  at  which  flame  blow-off  occurs;  d  Is  a 
determining  dimension;  un  is  the  normal  flame  velocity;  v  is  the  kine¬ 
matic  viscosity;  a  is  the  coefficient  of  thermal  diffuslvity;  K  Is  a 
constant. 
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HEAT  EXCHANGE  IN  THE  HORIZONTAL  SURFACES  OF 
HEAT  INSULATING  STRUCTURES  SLOTTED 
WITH  METAL  IN  DEEP  COOLING 
SYSTEMS 

N.  V.  Krylov  and  K.  I.  Strakhovich 


Heat  exchange  in  the  horizontal  surfaces 
of  heat  insulating  structures  slotted  with 
metallic  structures  in  systems  of  deep 
cooling  is  experimentally  investigated.  The 
values  of  Coefficient  of  heat  exchange  obtained 
are  compared  with  its  calculated  values  derived 
by  various  means. 


In  number  of  problems  of  contemporary  heat  technology  and 
cryogenics  it  is  necessary  to  determine  the  heat  exchange  on  the 
surface  of  different  types  of  the  heat  insulating  setups  slotted 
with  metal  structures.  The  existent  formulas  for  determining  the 
coefficient  of  heat  exchange  for  horizontal  and  vertical  walls  and 
cylindrical  surfaces  under  conditions  of  natural  convection  give 
values  for  this  coefficient  which  are  very  different  from  actuality, 
especially  at  low  temperatures  [1,  5]. 

This  is  explained  by  the  fact  that  a  purely  statistical 
approach  to  the  problem  of  determining  a  adopted  for  the  case 
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of  an  insulating  structure  slotted  with  metal  does  not  correspond 
to  the  actual  picture,  since  the  temperature  gradient  sharply 
increases  above  the  surface  of  the  metal  occupying  a  relatively 
small  area  in  the  whole  contact  surface  of  the  examined  insulating 
structure  with  air  or  cold  gas;  and  this  creates  a  corresponding 
temperature  pressure  in  the  contiguous  &ir  (gas)  layer  and  by  that 
very  fact  damages  the  pureily  statistical  conditions  of  heat  exchange 
via  the  given  layer  of  air  ( gas) . 

Hence  the  necessity  arises  of  investigating  the  distribution 
of  the  air  and  cold  gas  temperatures  above  the  insulating  surf dee 
and  above  the  surface  of  the  metal  in  order  to  Obtain  more  exaft 
values  for  coefficient  of  heat  exchange  a. 

The  experimental  arrangement  used  in  the  present  work  was  a 
metal  Dewar  vessel  of  diameter  d  *  250  and  of  height  h  =  600  mm 
(Fig.  l)  .  Into  this  vessel  was  put  a  laboratory  specimen  of  the 
insulating  structure  in  the  shape  of  a  cylinder  (b  *  200,  h  =  200  mm) 
which  was  filled  with  the  insulant  intersected  with  "heat  bridges" 
(metal)  .  These  heat  bridges  were  made  of  various  materials  and 
were  of  different  construction.  *j 

Depending  on  the  requirements  of  the  experiment  the  heat 
bridge  of  a  respected  material  either  was  supported  from  both  the 
hot  and  cold  side  on  steel  disks  ( diameter  *  200  mm  and  thickness 
fi  =  15  mm)  imitating  heat  barriers,  or  merely  cut  through  the 
insulation  without  heat  barriers.  Loose  insulation  (Mipor,  aerogel, 
and  slag  cotton) ,  was  utilized.  When  the  liquid  nitrogen  poured 
into  the  Dewar  vessel  reached  a  certain  level  and  when  the  temperatures 
throughout  the  vessel  and  the  insulation  surrounding  it  had  been 
distributed  in  a  constant  manner,  the  experimental  specimen  was 
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fixed  In  the  vessel  supported  on  stops  welded  to  the  internal  walls 
of  the  vessel.  The  vessel  meanwhile  was  fastened  In  a  hollow  steel 
cylinder  on  a  light  stand  and  the  clearance  between  its  external 
surface  and  the  inner  steel  cylinder  was  packed  with  insulation 
in  order  to  decrease  the  heat  currents  from  the  surrounding  medium 
to  the  liquid  nitrogen. 

In  order  to  determine  the  temperatures  in  the  mass  of  insulation 
along  the  heat  bridge  and  also  the  temperatures  of  the  air  bathing 
the  upper  surface  of  the  specimen, and  also  of  the  cold  gas  bathing 
the  lower  surface  of  the  specimen,  we  attached  copper- constantaji 
thermocouples  with  a  wire  diameter  d  *  0.2  mm  (Fig.  2). 

The  thermocouples  were  led  out  to  a  multipoint  switch  connected 
with  a  laboratory  potentiom  eter  for  accurate  measurements  of  the 
emf.  The  zero  point  was  secured  by,  placing  the-  warm  junctions  in 
a  glass  Dewar  vessel  filled  with  crushed  ice  and  dldtilled  wkt^r. 

The  control  indices  were  taken  off  until  a  sfcatjlbnery  regime 
was  reached.  Confirmation  of  this  was  the  constancy  of  the  emf 
values  in  the  respective  points.  ' 

From  processing  the  experimental  data  temperature  dlstrlbutld^is 
through  the  air  and  cold  gas  above  the  surface  and  insulation  and 
the  heat  bridges  were  obtained  as  a  function  of  the  material  from 
which  the  heat  bridges  were  made,  of  their  construction,  and  of 
the  methods  of  insulating  them  from  the  surrounding  medium  and 
heat  barriers  ( Fig.  2) .  The  adduced  data  show  that  in  a  definite 
area  of  the  air  at  a  height  of  about  20  to  40  milimeters  of  the 
heat  ba:.Tier  (insulation)  and  1  to  25  milimeters  above  the  surface 
of  the  heat  bridges,  whether  or  not  projecting  on  heat  barriers, 
there  is  observed  a  substantial  difference  In  the  temperature 


io- 
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distribution.  Above  the  indicated  region  the  same  air  temperature 
distribution  takes  place  for  all  the  examined  types  of  insulating 
mechanisms. 

Table  1  shows  the  relative  values  of  heat  coefficients  o 
calculated  by  Various  methods  for  the  heat  insulating  structure 
slotted  with  heat  bridges  constructed  as  shown  in  Fig.  2  and  made 
of  copper  and  lKhl8N9T  steel. 

For  heat  insulating  structures  traversedby  the  other  metals 
used  in  the  investigation,  the  values  of  the  coefficients  of  heat 
exchange  lie  within  the  limits  indicated  in  Table  1,  the  change  in 
a  occurring  in  direct  dependence  on  the  heat  conductivity  of  the 
metal  and  insulation. 

The  difference  in  the  coefficients  of  heat  exchange  o  from 
air  to  metal  calculated  by  the  indicated  methods  1b  explained 
by  the  fact  that  because  of  the  slight  thickness  of  the  boundary 
layer  it  is  very  difficult  to  measure  that  thickness  exactly;  the 
coefficients  of  heat  exchange  a,  however,  experimentally  obtained 
and  also  obtained  by  calculating  the  heat  currents  along  the 
thermo  elements,  differ  one  from  another  within  limits  permissible0 
in  our  calculation. 

The  increase  in  the  velocity  of  air  movement  on  the  warm  side 
of  the  insulation  structure  caused  by  a  substantial  temperature 
change  evidently  explains  the  nonconformity  of  the  experimental 
values  of  the  coefficients  of  heat  exchange  a  with  their  calculated 
values  obtained  from  extent  formulas. 

The  value  of  the  coefficient  heat  exchange  from  the  metal  to 
the  cold  gas  is  higher  in  comparison  with  its  value  from  air  to 
metal.  This  is  explained  by  the  formation  of  a  vapor  resulting 


from  the  evaporation  of  the  liquid  nitrogen  and  filled  with  the 
tiniest  droplets  of  liquid  whose  contact  with  the  warm  surface  of 
the  metal  leads  to  vigorous  vaporization  strengthening  the  heat 
exchange,  and  consequently  also  to  increasing  the  coefficient  of 
heat  exchange  a  from  metal  to  cold  gas. 

The  values  of  a  from  air  to  insulation  and  from  Insulation  to 
cold  gas,  based  on  experimental  data  and  determining  the  heat  currents 
along  the  elements  of  the  insulating  structure^  satisfactorily  agree, 
but  differ  greatly  from  the  values  calculated  from  the  extant  and 
empirical  formulas. 

This  considerable  deviation  is  connected  with  the  fact  that 
above  the  insulation  we  observe  air  stagnation,  confirmed  by 
approximate  calculation  of  the  velocity  of  air  motion  (0i22  to  0.26 
m/aec  above  the  metal  surface)  ,  but  the  ratio  of  the  velocity  of 
air  motion  above  the. metal  ut  and  above  the  insulation  u2  is  Ui/u*  ■ 
=  20-50.  It  is  obvious  that  at  such  small  velocities  their  In¬ 
considerable  change  leads  to  a  noticeable  change  in  coefficient  a; 
and  it  is  necessary  to  keep  this  in  mind  when  solving  practical 
problems.  ^ 


TABLE  1 


Relative  values  of  coefficients  of  heat  exchange  for 
heat  insulating  structures  traversed  with  heat  bridges. 
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8,15 

9.10 

9.0600 
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2.71 

6^72 

6.85 

0.0875 

air  to  metal  .  . 

5,600  ■ 

6.24 

8.12 

8.40 

9.20 

9,3000 

metal  to  cold  gas 

9,850 

2.84 

7,05 

7.10 

6.80 

2.4600 

lKhl8N9T 

Steel 

air  to  insulation 
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0,380 

5.54 

7.34 

7.65 

8.10 

8.3200 

insulation  to  col< 
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2.39 

5.92 

4.15 

m 
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Fig.  1.  Diagram  of  experimental  setup. 

I)  Steel  cylinder,  2)  specimen  under  in¬ 
vestigation,  5)  insulation,  4)  metal  Dewar 
vessel,  5)  Dewar  vessel  for  thermocouple  warm 
Junctions,  6)  multipoint  switch,  7) mirror 
galvanometer,  8)  metering  device,  9) power 
source,  10) battery  of  galvanic  cells* 

II)  normal  Weston  cell,  12)  potentiometer  Weston 
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63  115  M  179  211  213  275  T 


Fig.  2.  Change  In  temperature  of  air 
and  cold  gas  close  to  horizontal  surface 
of  Insulation  ( Mlpor)  and  of  heat  bridges* 
a) different  materials  in  heat  bridges  of 
the  same  construction:  air  side  ( 1-7) j 
cold  gas  side  ( 1  * -7  * ) J  horizontals 
( 1-7  and  1 1  —7 ' )  are  temperatures  of.  pro¬ 
jecting  surfaces  of  heat  bridges  of  copper 
1,  l’J;  aluminum  (2,  2’);  duraluminum 
Cu]  (3,  V) ;  brass  [7#Cu,  3(#Zn] 
k,  4*)i  carbon  steel  [0.01#  Cl  (5,  5 *)  i 
stainless  steel  [lKhl8N9T]  (6,  6')j 
textolite,  goethlnax(7/  7')  • 

(  b  and  c  continued  on  next  page) . 
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the  temperature  on  the  projections  of  which  Is 
shown  by  the  corresponding  numbers  on  the 
horizontals,  c)  Different  Insulation  structures 
at  the  ends  of  heat  bridges  all  of  the  same  material: 
1-5)  types  of  insulating  structures,  the  temperature 
on  the  projections  of  which  is  shown  by  the  corres¬ 
ponding  numbers  on  the  horizontals. 
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